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We study the routing of quantum information in qubit chains. This task is achieved by suitably 
chosen time-dependent local fields acting on the qubits. Employing the physics of coherent destruc- 
tion of tunneling, we demonstrate that a driving-induced renormalization of the coupling between 
neighboring qubits provides the key for controlling the transduction of quantum information be- 
tween permanently coupled qubits. We employ this idea for building a quantum router. Moreover, 
we discuss the experimental implementation with Penning traps, and study the robustness of our 
protocol under realistic experimental conditions, such as fabrication uncertainties and decoherence. 
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I. INTRODUCTION 

Entanglement, the inextricable quantum correlation 
between two physical objects, is a key resource for quan- 
tum informational tasks which promise pushing further 
our current computational limits This resource, 

which is produced through the interaction of two (or 
more) quantum objects, needs to be created, used and 
transferred under strict control such that a successful 
computation can be performed. Transporting entangle- 
ment may be considered as the analogue of electric cur- 
rents for classical computation, which need to be pro- 
duced, then used for computation inside a CPU, and fi- 
nally transferred to the output device. 

In this paper we will focus on the process of trans- 
ferring information once the computation has been per- 
formed. While long distance quantum communication 
can be done quite efficiently with photons [2J, commu- 
nication inside a quantum computer will probably use 
linearly arranged solid-state qubits [!, 0] • Lately, much 
attention has been paid to spin chains with nearest neigh- 
bor couplingj_which allow for perfect transfer of quantum 
information 0, 0, 0, HJ . The simplicity of such chains and 
the fact that they do not require special control make 
them robust for experimental implementations 9]. It is 
our belief that a quantum computer will profit from more 
complicated devices which allow not only linear propaga- 
tion, but also the possibility to stop the signal at a given 
node in a network, reverse its motion, etc., such that the 
information can be routed to the desired node with min- 
imal control. This could be particularly useful in quan- 
tum registers to store, retrieve, and transfer information 
contained in different qubits. 

The possibility of controlling the transfer in a spin 
chain is exciting. Such a chain could be used to inhibit 
access to parts of the register, or even to divide a signal 
into two different branches for distribution in a network 
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of nodes, provided that a particular propagation direc- 
tion can be chosen. At this point one could argue that 
for forbidding the propagation through a particular di- 
rection we "just need" to switch off the coupling in that 
direction [lfj. However, coupling and decoupling qubits 
is a formidable task [TJ [l2|, [ID, [lj] . Therefore we should 
think of an alternative for breaking the symmetry and, 
thus, allowing the transfer in one direction, without mod- 
ifying the inter-qubit couplings. 

The mathematics relating to coherent destruction of 
tunneling [l5|, |ji|[l7| or dynamical localization, i.e., the 
suppre ssion of ballistic transport by ac fields [13, [lg, Ixol . 
|20|.|21| suitably solves the objective in the high-frequency 
limit [22j |. Such driving renormalizes the tunnel matrix 
elements and has been investigated previousl y in the con- 
text of ac-driven mesoscopic transport (l8Ll2lL l23j |. In 
a recent theoretical work [2J|, it has been suggested to 
temporarily suppress tunneling between specific neigh- 
boring sites of an optical lattice, such that a particle will 
be transported into a direction of choice. Here we will 
demonstrate how to use such a dynamical renormaliza- 
tion for transmitting quantum information. 

In Sec.|TT]we introduce the concept of a quantum router 
and present in Sec. lIIII an according chain model together 
with brief explanation of the coherent destruction of tun- 
neling on which our quantum router relies. Section IIVI 
is devoted to the development of an according protocol 
for entanglement distribution and state transfer. A pro- 
posal for the implementation in Penning traps together 
with an analysis of the robustness of the method under 
experimental conditions is presented in Sec. El 



II. A ROUTER FOR QUANTUM 
INFORMATION 



Quantum information processing protocols usually de- 
pend on having full control of the qubits and their mutual 
interactions. This control relies on external devices such 
a lasers or signal generators, which limits the number of 
qubits by space considerations. One possible solution of 
this problem is splitting the processor into small blocks 
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FIG. 1: (Color online) Hub configuration. One qubit of an 
entangled pair is connected to a line with perfect transfer. 
Then the entanglement propagates and is later steered to one 
of the possible outputs, i.e., distributed to the desired trans- 
mission line. This enables quantum communication inside the 
quantum computer by both direct state transfer and by tele- 
portation once the entanglement of distant parties has been 
established. 



which, however, need to be connected [25|. Increasing 
the complexity of the network requires the selective cou- 
pling and decoupling of the blocks. We term such a cou- 
pler/decoupler quantum router. 

Let us assume for the moment that we have complete 
control of the coupling between every pair of qubits in an 
arbitrary array, despite that fact that this requires quite 
sophisticated engineering, at least much more than the 
perfect, unidirectional entanglement transfer described 
in fiefs. [1, H, 0|- Thus such a specialized piece of archi- 
tecture is probably useful only for very specific tasks at 
particular points during computation, transfer, and stor- 
age. One example is the hub sketched in Fig. [U let it be 
realized by a passive or an active element. When entan- 
glement arrives at the junction, it should be routed to 
an output of choice. The output channel may either be 
selected by the router itself or depend on an additional 
classical signal containing that information. Moreover, 
the router may select the output direction by the time at 
which the signal arrives at the hub. This means that the 
hub effectively opens the channels in turns and delivers 
the signal accordingly. 

A main obstacle for the experimental realization of a 
router is that it might be difficult if not impossible to sim- 
ply turn off the interaction between neighboring qubits. 
Still, there exists an indirect way for this task that is 
based on the physics of "coherent destruction of tunnel- 
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FIG. 2: (Color online) T-shaped configuration. Alice sends 
the state of one of her qubits without dispersion along the 
horizontal chain. When the state arrives at the node, an ac 
field with a ratchet-like profile transfers the signal to Bob or 
to Charlie. 



ing" . Let us consider a chain of qubits with given qubit- 
qubit interaction determined by the experimental setup. 
But nevertheless there exists the possibility of controlling 
the energy splitting of each qubit, which is typically less 
demanding than controlling the coupling. Then, as we 
will demonstrate below, an ac field with proper amplitude 
and frequency acting on the qubit will effectively renor- 
malize the qubit-qubit interaction. It is even possible to 
suppress this interaction almost entirely. After explain- 
ing the underlying mechanism, we apply the idea to the 
T-shaped setup sketched in Fig. [2l where Alice sends a 
quantum state through the horizontal quantum channel. 
When the state arrives at the node, it shall be directed 
either to Bob or to Charlie. We demonstrate below that 
this can be achieved by applying to the qubit splittings 
of the vertical chain an ac control with a ratchet-like am- 
plitude distribution. 



III. COHERENT DESTRUCTION OF 
TUNNELING IN CHAINS OF QUBITS 

We consider entanglement transfer in the vertical qubit 
chain of Fig. [^modelled by the Hamiltonian 

H =\Y. M*K + J £ to + < CT rt+l) (1) 

n n 

with o~n — \{o~n i icr*() and cr", a — x,y,z the usual 
Pauli matrices for spin n, while h = 1. This Hamil- 
tonian is known as the isotropic XY model with time 
dependent local fields h n (t). The qubit-qubit coupling 
is a standard interaction, realized either by direct engi- 
neering of the chain or as a rotating-wave approximation 
(RWA) where terms of the type cr+er+ +1 and cr~a~ +1 
are neglected [2(1 H3, UK]. By means of the Wigner- 
Jordan transformation, Hamiltonian {1} is equivalent to 
the one-dimensional spinless fermion model described by 
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the Hamiltonian 



H = \ Yl h n(t)cUn + (4cn+l + 4+1^) ■ ( 2 ) 



Furthermore, we will see that the mechanism explained 
below can be observed as well in a bosonic chain of har- 
monic oscillators for which the Hamiltonian reads 

1 ^ 

H = + muj l(t)x 2 n + X22x n x n+1 , (3) 



provided a RWA in the coupling between the oscillators 
can be applied, i.e., neglecting the counter- rotating terms 
a n a n+ i and o+o^ +1 . 

We assume that the onsite qubit splitting possesses a 
time dependence stemming from a classical field of the 
form 



h n (t) = Ab n cos(ojt) 



(4) 



In the interaction picture defined by X = [/q (i)X[/o(i), 
where Uo(t) = exp[— i(p(t)a z ] and ip(t) = h n (t)/2, the 
Hamiltonian (fTJ) reads 

H(t) = J^e^Wa+a- +1 + e~^a~a+ +1 (5) 

with the phase 

A 

(6) 



Vn(t) = — (b„ ~ & n +i)sin(wt). 



In the high-frequency limit lj J, the interaction- 
picture Hamiltonian can be replaced by its time- 
average. Thus, we obtain an effective static XY chain 
Hamiltonian for which the coupling between qubits n and 
n + 1 is renormalized according to 



J — > Jcff.ri — J Jo [A 



b n — b„+i 



(7) 



where Jo denotes the zeroth-order Bessel function of the 
first kind. Thus, the coupling between sites n and n + 1 
can be tuned to zero by choosing driving parameters for 
which A(b n — b n -\-i)/u> is a root of Jo- This renormaliza- 
tion has been used to explain the phe nomenon of coherent 
destruction of tunneling 0, [ltjl Il7l. l29j . Here we use it 
to steer quantum information by temporarily suppressing 
the interaction between two particular spins. 



IV. PROTOCOLS AND CONTROL 
MECHANISM 

Quantum routers may be useful for various purposes. 
We focus in the following on two of them, namely en- 
tanglement distribution and quantum state transfer. We 
develop for these aims protocols that rely on the CDT 
mechanism introduced in the previous section. It will 
turn out that both protocols are closely related. 



A. Entanglement routing 

An essential resource for quantum communication is a 
pair of distant entangled qubits, one at each end of the 
communication channel [301 ] . Such pairs can be obtained 
by creating the entangled pair locally and subsequently 
transfer to each end one partner. Recently, Creffield sug- 
gested to use the CDT mechanism described above for 
the controlled coherent propagation in a lattice [24[ ■ Dif- 
ferent to this entanglement transfer is the previously sug- 
gested non-local creation of entanglement in a spin chain 

The topic of this subsection, by contrast, is to create an 
entangled pair at Alice's place and to route one partner 
to either Bob or Charlie, see Fig. [5J A possible protocol 
for this task is the following: Let Alice have a two-qubit 
entangled state |V>)a that she wants to share with Bob. 
She will then attach one of the qubits to a chain with 
isotropic nearest neighbor interaction such that the qubit 
propagates towards the node with the vertical chain. If 
she for example owns the state \i/j(0)) = (|Ql) + |10))/\/2, 
the horizontal chain has the initial state 
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V2 



(|F M |00...00) + |0) j4 |10...00», (8) 



while the desired final state is 



1 



|V>)out = -^=(|1)a|00...00) + |0) a |00...01)). (9) 

Such perfect entanglement transfer can be achieved by 
specific static couplings [1, H, The question is now 
whether time-dependent fields allow one to route the 
qubit in a controlled manner from the node of the T- 
junction to either Bob or Charlie. 

When the travelling qubit arrives at the node of the 
T-junction (see Fig. [2]), coherent destruction of tunnel- 
ing must become active such that propagation to, say, 
Charlie's branch is suppressed. Our goal is now to find 
an ac field that such that the entanglement is perfectly 
transmitted to Bob. Moreover, the protocol should be 
flexible enough to allow routing to Charlie as well. 

We consider now a qubit chain that consists of blocks 
with four qubits with the ratchet-like energy splittings 

MM 



b, n = An' 

& + Ai, n = 4n' + l 

6 + Ai+A 2 , n = 4n' + 2 

b + A 2 , n = An' + 3 



(10) 



for integer n' , as is sketched in Fig. [31(a). It is constructed 
such that the energy differences between two neighbor- 
ing qubits are given by the sequence Ai, A2, — Ai, — A2, 
i.e., their absolute values alternate between Ai and A2. 
CDT can now be employed for temporarily suppressing 
tunneling between qubits with the one or the other value 
of the splitting [note that Jq(—x) = Jq(x), such that the 



4 



a) 



b+A+ A 2 



b+ A. 



b+ A, 



b) 



■3 




t[l/J] 



12 f 16 

Ti + r 2 



FIG. 3: (a) Spatial variation of the local qubit energy split- 
tings along the vertical chain. It is chosen such that the abso- 
lute value of differences between two neighboring sites alter- 
nates between Ai and A2. (b) ac field with alternating driving 
amplitude. The amplitudes [see Eq. are such that tun- 

neling between qubits with either energy difference Ai or A2 
is suppressed. 



CDT condition is not sensitive to the sign of Ai^]- A 
field that alternates between these two possibilities can 
induce directed transport [24J. We here use this idea to 
route entanglement. To be specific, we act on qubit n 
the driving field 



can be mapped to the tight-binding model 
H = J2 h n (t)\n){n\ + ^ £ (\n){n + l\ + \n+ l)(n\) , 

n n 

(12) 

with the (N + l)-dimensional state vector 



N 



(13) 



where |n) = |0)^|00 . . . 1„ . . . 00) denotes the state in 
which the nth qubit of the chain is the only one in state 
1. In order to achieve a compact notation, we have in- 
troduced the state |0) = |1)a|0 ... 0) for which all qubits 
are in the ground state, while Alice' state is excited. 

The main requirement for our protocol is that Alice' 
qubit remains entangled with the qubit that propagates 
to Bob or to Charlie, respectively. We measure this prop- 
erty by the concurrence C = max{Ai — A2 — A3 — A4, 0} 
|3fj |. The A's are the ordered square roots of the eigen- 
values of the matrix p(<rj a^)p*(u\ j <g> a^) with p be- 
ing the reduced density matrix of the considered qubit 
pair. In our case, the relevant reduced density ma- 
trix is that of Alice' qubit and qubit n of the chain. 
Tracing out all other qubits, we obtain in the basis 
{\0a0u}, \0a1u), |1a0«), I 1a In)} the expression 
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(14) 



for which the concurrence reads 



, , s , , , s fl/Ai, 0<t<Ti 

Kit) = Zo^bn COs(ujt) X ( 

I I/A2, J-l S t < 1\ + ±2 

(11) 

and periodically continued, see Fig. [3Jb), where £0 = 
2.4048 ... is the smallest positive root of the Bessel func- 
tion Jo. For this field, the CDT condition J c s — is ful- 
filled the transitions An' «-> An' + 1 and An' + 2 <-> An' + 3 
during the time interval [0, T\], while the other transi- 
tions are still allowed. For them the effective coupling 
energy is J e s,x — J</o(£oAi/A 2 ). The length of the time 
interval, 7\ — 7r/|J e fr,i|, corresponds to a half tunnel cy- 
cle between qubits An' and An' + 1, such that these qubits 
exchange their state. During time interval [Ti, T± + T2], 
we obtain that the coupling between spins with energy 
difference A2 is suppressed, while the other coupling is 
renormalized according to J c a^ = JJo(£oA2/Ai). Thus 
the second stage has to last for T2 = tt/| Jeff,2 1 - 

Next, we test numerically the quality of the entangle- 
ment transfer. Since the isotropic XY model conserves 
the number of excited spins, it is sufficient to consider 
only the subspace with at most one qubit in state 1, 
while all other qubits are in state 0. Then the chain 



Ca,u = 2|c c„|. (15) 

If Alice transfers her qubit directly to the node, i.e., to 
the central qubit ra no d e of the vertical chain, we can start 
our numerical calculation with the initial state 

|V(0)) = i(|0) + K odc )). (16) 

Figure |4] shows the corresponding time evolution of the 
concurrence Ca,u- It demonstrates that for a proper on- 
set of the driving field (fTTj) , the entanglement propagates 
solely to Bob who finally receives a state that is perfectly 
entangled with Alice' state. For a driving field starting 
with the larger of the two amplitudes, the entanglement 
propagates equally perfectly to Charlie, see Fig.O Thus, 
one can route the entanglement to a particular end of 
the vertical chain by choosing a proper initial time of the 
driving field. 

It is also possible that both Bob and Charlie receive a 
state that is entangled with Alice' state, as can be appre- 
ciated in Fig. [5] There the travelling qubit arrives with 
a delay of t = Ti/2, such that the routing towards Bob 
is only half completed at the time the upward channel is 
closed and routing to Charlie begins. 
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FIG. 4: (a) Time evolution of the concurrence between Alice' 
qubit and qubit n of the vertical chain with length N = 13 for 
the parameters uj = 10J, b = J, Ai = 3J, A2 = 1.5J. Alice' 
qubit is initially entangled with the qubit at n n ode = 7. The 
onset of the driving field (|11[) is such that initially transfer 
to Charlie is suppressed, (b) Same as panel (a) but for an 
arrival time t = Ti at the node, such that the entanglement 
propagates to Charlie. 



Our intuition for the dynamics of the quantum router is 
based on the derivation of an effective time-independent 
Hamiltonian derived within high-frequency approxima- 
tion. In the following we numerically test the results of 
our high-frequency approximation for finite driving fre- 
quencies. In Fig. [6] we plot the final concurrence Ca,b of 
Alice' qubit and the qubit that Bob receives, for different 
driving frequencies. We find that the protocol works, pro- 
vided that the driving frequency is sufficiently large. For 
chains of length up to N ~ 100, the final entanglement 
is quite large, e.g., Ca,b > 0.97 for w > 30J. For driving 
frequencies u> < 20 J, the concurrence already assumes a 
significantly lower value. 
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FIG. 5: Same as Fig. |4]for the arrival time t = Ti/2 at the 
node, such that the travelling qubit is split. 
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FIG. 6: Length dependence of the final concurrence for chains 
of various frequencies. The differences of the energy splittings 
are Ai = 100J and A2 = 59.02J, corresponding to the ratio 
A2/A1 = 0.5902 which provides the optimum transfer veloc- 
ity. 



In Fig. O we used a driving field with the amplitude 
ratio A1/A2 = 0.5902 for which the entanglement trans- 
fer is as fast as possible for a given value of J. Then 
the operation time and, thus, the total coherence loss is 
kept to a minimum, see Sec. IV A 21 In order to obtain 
the optimized amplitude ratio, one has to consider that 
during a full cycle of the driving field, the entanglement 
is transported over two sites. This takes the time 



Ti+T 2 = 



2|J (Ai/A 2 )| 2|J (A 2 /A 1 )|- 



(17) 



which possesses a minimum for the mentioned ratio and 
for the inverse ratio as well. 
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B. Quantum state transfer 

The spin chain can also be used for quantum state 
transfer, such that a classical communication channel be- 
comes dispensable @. For investigating state transfer, 
we assume that Alice owns the state 



\i/>) A = a\O)+0\l) 



(18) 



which she likes to transfer to Bob. In doing so, she ini- 
tializes the chain with the first qubit in the target state, 
while the other qubits of the horizontal chain are set to 
the ground state. For an excitation-preserving spin chain 
such as the one defined by the Hamiltonian (fTJ) , it can be 
shown that the unitary evolution of the chain provides 
Bob with the state Q 



Pb 



l-\f3\ 2 \c B \ a(3*c B 
a*(3c B \p\ 2 \c B 



(19) 



written in the basis {|0), |1)}. Perfect state transfer cor- 
responds to |cb| = 1. For our driven spin chain with 
different local qubit splittings and ac fields, each transfer 
from one qubit to the next contributes a particular phase 
shift. Thus, Bob will receive a state with cb = |cs|e l61 . 
In order to restore the full state (fT5)) , one has to know 
the total phase shift 6, which requires full knowledge of 
the transfer process, i.e., one has to calibrate the chain 
before use. Subsequently, a fixed unique phase shift has 
to be applied to every signal coming from the hub. Such 
hub together with a data sheet that contains the phase 
shifts could represent a useful piece of hardware. 

Finally, the quality of the state transfer can be quan- 
tified with the fidelity 



F = \{^ A \p B \^ A }\, 
which in our case reads 



(20) 



F = 1 - (1 - C A ,i 



[1-C a , b (1-2|/?| 2 )] , (21) 



after the phase has been rotated back. This means that 
the fidelity directly relates to the concurrence, C a ,b in 
Eq. (fT"5|) . of the entanglement router. 



V. PHYSICAL IMPLEMENTATIONS 

A highly accurate and stable confinement of ions in 
traps has been achieved in the last decade, together with 
the possibility of working in the ground state of the ion 
motion [36L \3l\ . This enables the realization of quantum 
chains that consist of harmonic oscillators with a nearest- 
neighbor coupling. The physical origin of the coupling 
can be Coulomb interaction between the confined parti- 
cles or a capacitive couplin g of the ions residing in neigh- 
boring wells; see e.g. Refs. [26l l27l] . For a comprehensive 
review of according control and readout techniques, see 
Ref. S3. 



The Hamiltonian for the excitations in Eq. ([2]), with 
a general onsite Zeeman field h n (t) can also be imple- 
mented by a chain of harmonic oscillators. The par- 
ticular form of the interaction is typically obtained by 
coupling the spatial coordinates of neighboring oscilla- 
tors [see Eq. (|3j)]. as is the case for Coulomb interac- 
tion between charged particles in dipole approximation. 
Note however, that the resulting Hamiltonian contains 
terms that create and annihilate phonon pairs. These 
terms vanish within rotating-wave approximation, which 
means that a mapping to our spin-chain model ([lj is 
possible only if this approximation holds, i.e., if the cou- 
pling is much weaker than the time average of an applied 
"Zeeman" field, J -€i ho, where here h n (t) = ho + b n {t). 
Notice that the presence of the constant "Zeeman" field 
ho merely rescales the energies and will not affect the 
router otherwise. Experimentally, this field can be re- 
alized by an ac voltage with non-zero mean, such that 
h n {t) — 2htu n (t) with the frequency w n (t) proportional 
to the square root of the applied voltage. 

As a particular example, we consider singly trapped 
ions in planar Penning traps, in which neighboring ions 
couple capacitively through a conducting wire. Neglect- 
ing resonant absorption in the wire, the coupling for a 
typical trap size of 1mm will be roughly 100 times lower 
than the onsite potential [27| . For even smaller traps 
of size 0.1mm, these quantities still differ by a factor of 
10. Thus, engineering the coupling strength essentially 
means choosing the correct trap size. For planar traps, 
the ion can be moved further away from the electrodes, 
which reduces the coupling even more. In a linear ar- 
ray of trapped ions such as the one of Ref. [26[, the 
atoms experience Coulomb interaction. In both cases, 
the ratio between interaction strength and onsite poten- 
tial hui is proportional to q 2 /(ireomuj 2 d 3 ). The difference 
in strength between both approaches is a form factor, 
given by the fact that the wire-mediated coupling has a 
given efficiency when transmitting Coulomb interaction. 
In the case of direct Coulomb interaction, the coupling 
is slightly stronger, but can be reduced by pushing the 
ions far apart from each other, using that fact that J is 
reciprocal to the typical distance between the ions. 

It must be noted here that for Penning traps, all pa- 
rameters basically depend on two scales, namely the trap 
size and the voltage. The motional frequency is propor- 
tional to W jd, with V the applied voltage and d the 
size of the trap. Manipulating both scales allows one to 
achieve the desired parameter values. 



A. Robustness under realistic conditions 

In an experiment, one usually faces additional difficul- 
ties not captured by a Hamiltonian with perfectly stable 
parameters. Two such difficulties come to mind, namely 
imperfections in the fabrication process and the unavoid- 
able influence of an environment that causes dissipation 
and decoherence. 



1. Fabrication uncertainties 

An essential ingredient to our entanglement distri- 
bution protocol is the action of the onsite ac fields 
(t;oui/Ai) smcot during times Ti = ir/J e s^ 
with J effi i( 2 ) = JJo(CoAi(2)/A 2( i)), [see Eq. (TTJ) and 
Fig. [3] . This implies two critical experimental require- 
ments. First, the driving amplitudes have to match with 
good precision the first zero of the Bessel function and, 
second, the field amplitudes have to be switched after a 
time T12. 

Trap diameters in the sub-millimeter range, corre- 
spond to motional frequencies of the order GHz, with 
a J ~ 100MHz, such that the rotating wave approxima- 
tion applies. Moreover, a driving frequency of the order 
0.1 GHZ is available with current technology, while com- 
mercially available function generators at those frequen- 
cies can have an accuracy of several 10 9 samples per sec- 
ond. Thus the switching times can be adjusted with pre- 
cision tT — 10 _3 /J. The uncertainty in J mostly comes 
from measuring the distance between trapped ions (typ- 
ically by fluorescence). Assuming an accuracy of 10%, 
we obtain the error ej = 10~ 7 J, which corresponds to 
£t = 10 -7 / J. Thus, the relevant restriction is the men- 
tioned uncertainty in the function generator, such that 
€t = 10~ 3 /J. Concerning the amplitudes, we assume 
a relative error in the applied confining voltages of the 
order of 10 -6 , which yields an amplitude uncertainty of 
the order of 10~ 7 J. 

For a numerical simulation of these errors, we choose 
for each coupling matrix element Ji a random value from 
the interval [J — ej, J + e,j] with equal distribution. The 
interval width 2ej is determined by the experimental un- 
certainties discussed above. Accordingly, we select a field 
amplitude from the interval [b n — ef,, b„ + e&]. In the fol- 
lowing, we denote ensemble averages of such realizations 
by an overbar. 

We already emphasized in Sec. IIIII that the protocol 
for entanglement transfer is closely related to the one for 
state transfer. The main difference is that the latter re- 
quires knowledge of the extra phase 9. It has to be gauged 
or determined from theoretical considerations, such that 
Bob is able to compensate this phase by a local transfor- 
mation. When the system parameters acquire a random 
component, the extra phase becomes random as well. 

In principle, the additional random phase can be com- 
pensated by local transformations. Thus, it cannot influ- 
ence the entanglement. Numerically, we have tested that 
the experimental uncertainties mainly affect this phase, 
making the concurrence fairly independent of this type 
of errors. Therefore we focus on the impact of the fab- 
rication uncertainties on the state transfer. Assuming 
that we can gauge our apparatus, we know that the ac- 
tual value of 9 = 9 + Ag deviates from the average by 
Ag. Notice that 9 can be obtained just by averaging over 
many realizations, as we discuss below. Then the fidelity 
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FIG. 7: Fidelity distribution for the initial state \i/j)a = 
(|0) + \1))/V2 in a chain with length TV = 40 and the driving 
frequencies (a) u> = 30 and (b) uo = 40. The qubits splittings 
are Ai = 100J and A2 = 59.02J, while the parameter uncer- 
tainties read ej = et, = 10 -7 and er = 10 -3 . The data have 
been obtained from 500 realizations. The resulting standard 
deviations of the concurrence are 0.01 and 0.1, respectively. 



(l2Tj)) becomes [cf. Eq. (1ST]) ] 

F=l-\f3\ 2 (l-2C A , B cosA e + CX B ) 
+ 2\l3\ 4 C A ^ B (C A , B -cosA e ) . 

Thus, all the possible errors are captured by the phase 
deviation 9 = 9 + Ag. 

Figure [7] shows the frequentness of the accordingly re- 
duced values for the fidelity. We have computed the fi- 
delity (|20[) for 500 realizations of the spin chain, and 
plotted a histogram with bin size 10~ 3 . Obviously, the 
driving frequency u> = 30 J yields in most cases a fidelity 
F > 0.98 [panel (a)]. For a larger driving frequency 
[panel (b)], the fidelity distribution is still peaked at a 
value larger than 0.99. The standard deviation, however, 
is quite large, such that significantly smaller fidelities be- 
come rather likely. Thus, fabrication uncertainties im- 
pose that the driving frequency should not exceed 30J. 
On the other hand, we have seen in Sec. IIV AI that this 
value in fact also represents a lower bound for the validity 
of the rotating-wave approximation within which the ef- 
fective coupling strength ([7]) has been derived. Thus, we 
conclude that the results are best for a driving frequency 
of the order lo = 30 J. 
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2. Dissipation and decoherence 

Still we have to consider the main obstacle for any 
quantum information task, which is dissipation and deco- 
herence caused by weak but unavoidable coupling of the 
processor to a typically huge number of uncontrollable 
degrees of freedom. Irrespective of its physical nature, 
the environment is usually modeled as a bath of harmonic 
oscillators, where each oscillator couples via its position 
operator to a system coordinate [Tt], [H, HtJ H(| ■ Here we 
assume that each qubit n undergoes pure and indepen- 
dent dephasing, i.e., that it its coordinate couples to 
a separate bath. 

For weak dissipation, one can eliminate the bath within 
second-order perturbation theory and derive a Bloch- 
Redfield master equation for the reduced density oper- 
ator p of the qubits Considering, only phase noise 
on the traps the equation takes a Lindblad form [42j , 

N 

p = - 1 j[h,p\ + \y.«p<-p)> ( 23 ) 

n=l 

where we have assumed that the effective decoherence 
rate 7 is the same for all qubits. Remarkably, such a 
master equation describes, apart from a transient, the 
decoherence process quite accurately [43|. For qubit ar- 
rays implemented with calcium ions confined in a Paul 
trap, the dephasing rate is 7 ~ 1 kHz, while the qubit 
splitting is several MHz. The heating rate, which would 
include relaxation terms in the master equation, is 190 
times smaller [44j. This justifies the phase noise model 
used here. We will use this numbers for a rather conserva- 
tive estimate and assume 7 ~ 10~ 4 J, where J ~ 10 MHz. 

For qubit chains under the influence of phase noise, the 
entanglement decay is determined by the coherence loss, 
and one finds that the concurrence decays exp onentially 
with a rate 27, i.e., Ca.bW ^ exp(— 27^) [45(]. This can 
be verified easily by considering the action of the dissi- 
pative kernel in the master equation on the off-diagonal 
density matrix elements p nn > . The entanglement propa- 
gation from the middle of the chain to Bob's end takes 
the time tab = (Ti + T 2 )N/4: [Recall that our ratchet 
mechanism transports signals during a time T\ + Ti over 
two sites; see derivation of Eq. (JT7J)] . Thus, we expect 
that decoherence reduces the finally achieved entangle- 
ment according to 

Cam =C AB (7=0)e- 2 ^ B . (24) 

We numerically computed the time evolution described 
by the master equation (|23p . The results shown in 
Fig.[5]confirm our expectations for the exponential decay. 
Moreover, they demonstrate that a decoherence rate of 
7 = 10~ 4 J reduces the concurrence in a chain of length 
N = 60 from Ca,b{j=0) = 0.98 in the absence of deco- 
herence to a value Ca,b = 0.95. This is still sufficiently 
large for most possible applications. For a larger deco- 
herence rate, 7 = 10~ 3 J , however, the concurrence drops 
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FIG. 8: Concurrence decay as a function of the chain length 
for two different decoherence rates 7 (symbols) and u — 30. 
The solid lines mark the estimate (|24p . All other parameters 
are as in Fig. [7] 

significantly. For example in a chain of length N — 60, we 
have the already quite small value Ca,b = 0.65. Thus, we 
can conclude that for chains of this length, transmission 
of sufficient entanglement is only possible for decoherence 
rates that fulfill 7 < 10~ 4 J. For stronger decoherence, 
a hub with shorter chains is nevertheless possible, in the 
extreme limit even one with branches that consist of just 
one qubit. 

We close our analysis of fabrication uncertainties and 
decoherence by addressing also the fidelity of the state 
transfer protocol. In the presence of decoherence, it is no 
longer possible to establish an exact relation between the 
fidelity and the concurrence [cf. Eq. (|2"Tj) ]. The reason 
for this is that state transfer is limited by the dephas- 
ing of individual qubits, while the entanglement decay 
is influenced by the coherence between different qubits. 
The state transfer fidelity to Bob can still be expressed 
in terms of Bob's reduced density matrix (|19p and now 
reads 

F = 1 - pu + \/3\ 2 (2p u - 1) + 2|pi | Re(a/3e iA *), (25) 

which implies that again the average phase shift 6 plays 
a role; cf. Eq. ([22)). 

We find numerically, that the final population of Bob's 
qubit, which is described by the density matrix element 
Pbb, is practically independent of both the decoherence 
strength 7 and the fabrication uncertainties. Without 
these influences, the state transfer protocol works al- 
most perfectly, such that we can approximate this den- 
sity matrix element by pbb = \P\ 2 ■ Moreover, we al- 
ready argued that off-diagonal matrix elements decay as 
Pbo — ol(5* exp(— 77-^45) and, thus, the finally achieved 
fidelity can be well aproximated by 

F = 1 — 2(|/3| 2 - |/3| 4 )[1 - cos(A e )e-^ B ] . (26) 

Thus, within some reasonable approximations, we again 
find for the fidelity a closed expression that is a function 
of the initial state. In fact, for 7 = and Cab (7 = 
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FIG. 9: Fidelity spread for a chain of length N = 40 with 
the initial state \iP}a = (|0) + jl))/V2- The driving frequency 
is 10 — 30 and dissipation strengths 7 = 10~ 3 J (a) and 7 = 
10" 4 J (b). The histo gram is obtained from 500 realizations 



yielding a standard deviation of 0.01 in both of them, 
other parameters are as in Figs. [7] [5] 



All 



0) = 1, expression (|2"rj)) becomes identical with the fi- 
delity (|22]) for an ideal chain. This also demonstrates 
that the fidelity is significantly affected by both fabrica- 
tion uncertainties and decoherence, expressed in Eq. (|26p 
by the random phase shift Ag and the exponentially de- 
caying factor, respectively. 

It remains to numerically corroborate the quality of the 
approximate result ([26]) . Figure [9] depicts the numerical 
results for a chain of length N = 40. It reveals that the 
maximum fidelity in the histogram coincides with making 
A e = in ([H]), namely F = 0.99 for 7 = 10~ 4 J and 
F = 0.93 for 7 = 10~ 3 J. Further the data spread, due to 
the fabrication errors, is pretty the same that in the case 
without decoherece [cf. Fig. \7\. This underlines that 
state transfer with fidelity F > 0.97 is achievable with 
the experimental realization discussed above, provided 
the dissipation rate does not exceed the value 7 = 10~ 4 J. 



VI. DISCUSSION AND CONCLUSION 

We have proposed a quantum router that can be em- 
ployed for steering one qubit of an entangled pair from 
Alice to either Bob or Charlie. A second, related appli- 
cation is quantum state transfer. With the experimental 



realization with ion traps in mind, we have presented 
analytical and numerical results that underline the feasi- 
bility of our scheme. 

The underlying physics is coherent destruction of tun- 
neling between two neighboring qubits. This can be 
achieved by ac fields with a certain ratio between driving 
amplitude and frequency. Since the effective amplitude 
depends on the difference of the qubits' Zeeman ener- 
gies, a proper sequence of these energies temporarily sup- 
presses for each qubit the interaction with exactly one of 
its two neighbors. After half a tunnel cycle, the field am- 
plitude is switched, such that the suppression is active for 
the connection to the other neighbor. As a consequence, 
the qubit is transported into one particular direction. A 
main benefit of this mechanism is that it does not rely 
on the local control of interaction parameters, but only 
on the less demanding application of proper ac fields. 

In the T-shaped configuration under investigation, one 
partner of an entangled pair is first routed to the node 
that connects Alice with Bob and Charlie. When the 
qubit arrives at the node, the further propagation de- 
pends on which branch is temporarily closed at that 
very moment. In particular, we demonstrated that for 
a proper phase of the driving field, the qubit is reli- 
ably routed towards the one or the other direction. The 
achievable entanglement increases with the driving fre- 
quency. For an idealized chain of 100 qubits with realistic 
parameters, the corresponding concurrence may assume 
values of the order 0.98. 

In an experiment, however, the setup and the opera- 
tion will not be ideal. Therefore, we have investigated 
two limiting influences, namely fabrication imperfections 
of the desired splitting sequence and the influence of ex- 
ternal degrees of freedom which lead to decoherence. The 
former reduce the achievable fidelity, in particular for 
high driving frequencies, i.e., in the regime in which the 
idealized model works almost perfectly. The numerical 
simulation revealed the existence of an intermediate fre- 
quency regime, in which realistic fabrication errors do 
not significantly reduce the achievable fidelity. 

The second limitation leads to the ubiquitous deco- 
herence which generally represents the main obstacle for 
the implementation of quantum information schemes. It 
turned out that dephasing does not limit the entangle- 
ment or state transfer as long as its rate is roughly three 
orders of magnitude smaller than the typical qubit split- 
ting. 

In summary, we have demonstrated that the proposed 
quantum router should work under realistic conditions 
such as fabrication uncertainties and decoherence. The 
implementation is feasible not only with ion chains, but 
with practically all interacting qubit arrays in a T-shaped 
configuration. Thus our quantum router may represent 
an essential building block in future quantum information 
networks. 
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